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SUMMARY 


The  influence  of  a  small  hole  at  the  tip  of  a  crack  on  the  compliance 
of  a  cracked  specimen  is  examined.  It  is  shown  that  the  compliance  of  the 
specimen  with  a  hole  is  the  same  as  the  compliance  of  a  specimen  with 
an  infinitely  sharp  but  slightly  longer  crack,  the  increment  of  length  being 
proportional  to  the  radius  of  the  hole.  In  the  case  of  mode  III  straining 
(the  tearing  mode)  the  equivalent  sharp  crack  extends  beyond  the  edge  of 
the  hole  a  distance  equal  to  the  radius  of  the  hole.  The  results  are  pertinent 
to  the  experimental  measurement  of  stress  intensity  factors. 


SOMMAIRE 


On  etudie  l’influence  d’un  petit  trou  a  fond  d’une  fissure  sur  la 
raideur  d’une  4prouvette  fissuree.  On  demontre  que  la  raideur  de  l’eprou- 
vette  est  egale  a  celle  d’une  6prouvette  avec  une  fissure  infiniment  aigue  mais 
un  peu  plus  longue,  le  prolongement  etant  proportionnel  au  rayon  du  trou. 
Dans  le  cas  du  mode  III  d’ouverture  de  la  fissure  (mode  de  dechirement)  la 
distance  entre  la  pointe  de  la  fissure  aigue  Equivalent*  et  le  hord  du  trou  est 
egale  au  rayon  du  trou.  Les  resultats  ont  rapport  a  la  mesure  expert mentale 
des  facte urs  d’intensite  des  contraintes. 
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THE  EFFECTIVE  LENGTH  OF  A  CRACK  TERMINATED  BY  A  HOLE 


1.0  INTRODUCTION 

It  is  a  well-known  result  of  fracture  mechanics  that  the  stress  intensity  factor  for  a  cracked 
specimen  can  be  related  to  the  change  in  compliance  of  the  specimen  as  the  crack  extends.  This 
result  provides  the  basis  for  experimental  determinations  of  stress  intensity  factors.  The  determination 
requires  measurement  of  the  specimen’s  compliance  for  a  range  of  crack  lengths  and  the  calculation 
of  rate  of  change  of  compliance  from  the  data. 

When  stress  intensity  factors  are  determined  in  this  way,  it  is  a  common  experimental 
practice  to  terminate  the  crack  with  a  small  drilled  hole.  Use  of  a  hole  facilitates  measurement  of 
crack  lengths  by  providing  a  smooth  well-defined  edge  whose  location  is  easily  set  and  measured. 
There  is  the  further  advantage  that  the  hole  blunts  the  crack  tip,  thus  reducing  the  danger  that  the 
crack  will  extend  when  the  specimen  is  loaded  for  the  measurement  of  compliance. 

Use  of  a  hole,  however,  has  the  disadvantage  of  altering  the  geometry  of  the  crack  tip  and 
the  neighbouring  stress  field.  Since  the  stress  intensity  factor  is  a  theoretical  concept  based  on  the 
idealization  of  an  infinitely  sharp  crack,  it  is  not  immediately  clear  how  experimental  measurements, 
obtained  from  a  crack  terminated  in  a  hole,  should  be  interpreted.  One  may  reasonably  assume  that 
the  perturbation  introduced  by  the  hole  is  small  provided  that  the  diameter  of  the  hole  is  sufficiently 
small,  but  it  is  difficult  to  say  anything  more  precise  without  further  analysis. 

The  present  note  offers  an  analysis  of  the  change  in  compliance  caused  by  the  presence 
of  a  hole  at  the  end  of  a  crack.  It  is  found  that  the  compliance  of  the  specimen  with  a  hole  is  the 
same  as  the  compliance  of  a  specimen  with  an  infinitely  sharp,  but  somewhat  longer,  crack.  The 
increment  of  length  between  the  crack  ending  in  a  hole  and  the  equivalent  sharp  crack  is  found  to  be 
directly  proportional  to  the  radius  of  the  hole,  with  the  constant  of  proportionality  being  essentially 
the  same  fixed  quantity  for  all  configurations  of  specimen. 

Although  the  constant  of  proportionality  is  independent  of  specimen  configuration  it  is 
not  independent  of  mode  of  straining.  Different  constants  apply  to  the  three  basic  modes  of  straining, 
that  is  to  the  opening  mode,  the  sliding  mode  and  the  tearing  mode.  Determination  of  the  constants 
requires  an  analysis  of  the  stress  field  surrounding  the  hole.  The  analysis  is  carried  out  in  this  report 
for  the  tearing  mode  but  not,  because  of  mathematical  difficulties,  for  the  more  important  opening  and 
sliding  modes.  In  the  case  of  the  tearing  mode  it  is  found  that  the  equivalent  sharp  crack  extends  beyond 
the  edge  of  the  hole  a  distance  equal  to  the  radius  of  the  hole. 

The  present  results  complement  previous  results  of  Bowie  and  Neal  (Ref.  1)  for  shallow 
edge  notches,  and  supercede  a  proposal  of  Srawley,  Jones  and  Gross  (Ref.  2)  for  equivalent  crack 
length  which  is  now  seen  to  be  erroneous. 


2.0  ANALYSIS 

2.1  Comparison  of  Real  and  Ideal  Specimens 

Consider  two  cracked  specimens  A  and  B  of  general  shape,  as  shown  in  Figure  1.  Specimen 
A  is  the  ideal  specimen  in  which  the  crack  tip  is  infinitely  sharp,  while  specimen  B  in  the  real  specimen 
in  which  the  crack  is  terminated  by  a  small  hole.  Except  for  the  different  configuration  at  the  crack 
tip  the  two  specimens  are  of  identical  shape  and  both  are  loaded  in  the  same  way.  The  task  is  to 
relate  the  compliance  of  the  ideal  specimen,  from  which  the  crack-tip  stress  intensity  factor  can  be 
determined,  to  the  compliance  of  tile  real  specimen,  which  is  the  quantity  that  is  actually  measured 
experimentally.  Since  the  compliance  in  this  context  is  essentially  a  measure  of  the  strain  energy  in  a 
body,  the  task  is  equivalent  to  obtaining  a  relation  between  the  strain  energies  of  the  real  and  ideal 
specimens. 
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If  the  stress  field  Oy  which  exists  in  the  ideal  specimen  A  were  applied  to  the  real  specimen 
B,  it  would  satisfy  the  boundary  conditions  on  the  external  boundary  and  on  the  crack  but  would 
leave  residual  tractions  T*  on  the  surface  of  the  hole.  The  residual  tractions  T*,  which  are  self- 
equilibrating,  can  be  annulled  by  superposing  a  stress  field  a  y  which  leaves  tractions  -  T*  on  the 
surface  of  the  hole  and  which  satisfies  zero  boundary  conditions  elsewhere.  Thus  the  stress  field 
a\ j  in  the  real  specimen  B  is  the  superposition  of  fields  Oy  and  o'j ; 

ail  =  ay  +  a[s  (1) 

If  €jj,  e-j ,  €jj  are  the  strains  correcponding  to  stresses  Oy ,  Oy,  ojj  respectively,  the  strain 
energy  density  in  the  real  specimen  is 

^  tt  H  _  ^  ,  .  #. 

2°ijeU  =  2  + 


1  1  ,  , 

2  °U  eij  +  £  °'i  e‘i  +  a'i  e>i 


(2) 


where  we  have  used  the  fact  Oy  ey  =  a\.  ey  which  is  a  consequence  of  the  linear  relation  between 
stress  and  strain.  The  strain  energies  of  tne  specimens  are  obtained  by  integration  of  the  strain  energy 
densities  over  the  regions  occupied  by  the  specimens.  It  is  convenient  to  define  two  regions  Rj  and 
R2  relative  to  either  specimen,  where  R,  is  the  region  exterior  to  the  boundary  of  the  hole  and  R2 
is  the  region  within  the  boundary  of  the  hole.  Thus  the  ideal  specimen  A  occupies  the  combination  of 
regions  R{  and  R2  while  the  real  specimen  B  covers  region  R,  only.  The  strain  energies  UA  and  UB 
of  specimens  A  and  B  therefore  are 


UA  =  /  eadR  +  f~°aei)dR 

n  Z  D  “ 


Uu  = 


tt 

o  aM  e'i 


dR 


(3) 

(4) 


Now  integration  of  (2)  over  region  R,  and  substitution  from  (3)  and  (4)  yields 


The  first  integral  in  (5)  can  be  transformed  in  the  usual  way  into  an  integral  around  the  boundary 
of  region  R2  of  the  product  of  the  surface  tractions  and  displacements  associated  with  the  stress 
field  Oy .  Since  this  field  leaves  residual  tractions  Tf  on  region  Rj ,  then  it  leaves  equal  and  opposite 
tractions  -  T*  on  region  R2 .  Hence 

/  T°ij  «ij  dR  *  "  #4  T?ui  ds  <6> 

R2  *  Hz 

where  U|  is  the  displacement  field  associated  with  Oy  and  where  H  denotes  the  boundary  of  the  hole. 
The  third  integral  in  (5)  can  be  transformed  in  the  same  way  into  an  integral  around  the  boundary 
of  region  R| .  Since  the  tractions  associated  with  the  stress  field  Oy  are  zero  everywhere  except  on 
the  boundary  of  the  hole,  where  they  are  equal  to  -T*  we  have 


/  °!j  eij  dR  m  ~  i  T*Uj  ds 
R,  H 


(7) 
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It  follows  from  (5),  (6)  and  (7)  that 

UB  =  ua  +  /  ^  %  e»  dR  +  /  ~  o|j  ejj  dR  (8) 

R2  2  R,  2 

Formula  (8)  gives  the  relation  between  the  strain  energies  UA ,  and  UB ,  and  it  will  be  useful 
if  the  two  integrals  can  be  easily  evaluated.  Note,  however,  that  the  analysis  so  far  applies  to  all  three 
modes  of  straining. 

2.2  Dimensional  Form  of  the  Integrals 

Evaluation  of  the  first  integral  in  (8)  is  straightforward.  Since  the  integration  extends  only 
over  region  R2  which  covers  the  immediate  vicinity  of  the  crack-tip,  the  stress  field  Oy  in  R2  is  given 
with  adequate  accuracy  by  the  standard  expressions  for  stresses  near  a  crack  tip.  The  energy  density 
can  be  calculated  from  these  expressions  and  the  integral  evaluated.  Dimensionally,  the  result  will  be 
of  the  form 


/  \  ou  ey  dR  =  f,  K\  b/M  (9) 

R2  * 

for  a  specimen  of  unit  thickness,  where  KA  is  the  stress  intensity  factor  for  the  ideal  specimen  A,  6  is 
the  radius  of  the  hole,  f,  is  a  numerical  factor  and  M  is  the  appropriate  elastic  modulus.  For  mode  I 
or  mode  II  straining  M  is  Young’s  modulus  E,  while  for  mode  III  straining  M  is  the  shear  modulus  G. 

Evaluation  of  the  second  integral  in  (8)  depends  on  first  finding  the  stress  field  Oy  by 
means  of  a  stress  analysis.  It  is,  however,  possible  to  establish  a  priori  the  dimensional  form  of  the 
result.  Since  the  edge  of  the  hole  is  in  the  immediate  vicinity  of  the  crack  tip,  the  boundary  tractions 
-T*  which  determine  the  stress  field  can  be  obtained  with  adequate  accuracy  from  the  standard 
expressions  for  stresses  near  a  crack  tip.  The  tractions  are  therefore  proportionto  to  KA /y/b  and 
hence  the  stress  field  Oy  is  also  proportional  to  the  same  factor.  The  strain  energy  density  then  is 
proportional  to  K^/bM.  The  integral  of  strain  energy  density,  for  a  specimen  of  unit  thickness,  will 
be  proportional  to  /b M  times  the  square  of  some  characteristic  length.  Since  the  stress  field  Oy  is 
no  doubt  localized  around  the  hole,  the  hole  radius  b  is  an  appropriate  characteristic  length,  so  that 
the  second  integral  in  (8)  will  be  of  the  form 

f  \  Oy  ey  dR  =  f2  K\  b/M  (10) 

Rl  * 

where  f2  is  a  numerical  factor.  Because  of  the  localization  of  the  stress  field  Oy  around  the  hole,  the 
factor  f2  can  be  expected  to  be  largely  independent  of  the  overall  shape  of  the  specimen. 

The  relation  (8)  between  the  strain  energies  of  the  real  and  ideal  specimens  then  becomes 

UB  =  Ua  +  (*i  +f2>Ki  6/M  (11) 


2.3  Equivalent  Crack  Length 

The  stress  intensity  factor  for  the  ideal  specimen  is  related  to  the  strain  energy  release  rate 
by 


do  "  2M 


(12) 
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where  a  denotes  the  crack  length.  The  measured  stress  intensity  factor  KB ,  that  is  the  stress  intensity 
factor  determined  from  compliance  measurements  on  the  real  specimen,  is  in  effect  calculated  from 
the  analogous  expression 


dUB  _  k£ 
da  '  2M 

By  differentiating  (11)  with  respect  to  a,  and  using  (12)  and  (13)  we  find 

K|  =  K\  +  4(f,  +  f2 )  6  Ka  (dKA  /da) 


(13) 


(14) 


In  performing  the  differentiation  we  have  assumed,  as  will  be  verified,  that  f2  is  not  a  function  of  the 
crack  length  a.  This  assumption  is  reasonable  in  view  of  the  expected  concentration  of  the  stress  field 
oij  around  the  hole.  Taking  the  square  root  of  (14)  gives 

Kb  =  Ka  { 1  +  4  (f,  +  f2 )  b  (dKA  /do)/KA  )V2  (15) 

The  second  term  in  the  curly  bracket  is  no  doubt  a  small  correction  so  it  should  be  legitimate  to 
approximate  the  bracket  by  the  first  two  terms  of  its  binomial  expansion.  Thus 

Kb  =  KAjl+2(f,  +f2)ft(dKA/da)/KA[  (16) 


or 


Kb  =  Ka  +  2  (f,  +  f2 )  b  (dKA  Ida)  (17) 

which  is  the  formula  relating  the  exact  stress  intensity  factor  KA  to  the  experimentally  measured 
factor  KB . 

The  geometrical  interpretation  of  (17)  is  shown  in  Figure  2.  It  may  be  seen  that  if  KA  is 
the  stress  intensity  factor  for  a  crack  of  length  a,  then  KB  is  the  stress  intensity  factor  for  an  ideal 
crack  whose  length  is  greater  than  a  by  the  amount  2  (f  j  +  f2 )  b.  In  other  words,  terminating  a  crack 
by  a  small  hole  increases  its  effective  length  by  the  amount  2  (fj  +  f2 )  b. 

The  next  step  in  the  analysis  is  to  determine  the  values  of  the  factors  f ,  and  f2 . 

2.4  Factors  f(  and  f2  for  Mode  III  Straining 

The  preceding  analysis  applies  to  all  three  modes  of  crack  straining.  It  will  be  recalled  that 
the  three  basic  modes  of  straining  are  the  crack  opening,  also  known  as  mode  I;  the  sliding  mode  also 
known  as  mode  II;  and  the  tearing  mode,  or  mode  III.  From  now  on  the  analysis  will  deal  with  mode 
III  only.  Attention  is  confined  to  the  tearing  mode  solely  because  the  mathematics  involved  are  simple 
enough  to  cope  with.  The  mathematical  analysis  required  for  modes  I  and  II  is,  by  comparison,  far 
more  difficult  and  will  not  be  attempted  at  this  time. 

Up  to  now  the  precise  location  of  the  centre  of  the  hole  has  not  been  specified.  The 
preceding  analysis  is  valid  so  long  as  the  tip  of  the  crack  lies  somewhere  within  the  boundary  of  the 
hole.  It  is  now  necessary  to  fix  the  location  of  the  hole,  and  for  convenience  we  take  the  centre 
of  the  hole  to  lie  at  the  point  corresponding  to  the  tip  of  the  crack  in  the  ideal  specimen.  We  also 
introduce  a  system  to  coordinates  with  origin  at  the  centre  of  the  hole  and  oriented  so  that  the  crack 
is  aligned  with  the  negative  x-axis  as  shown  in  Figure  3. 

2.4.1  The  Factor  f, 

For  mode  III  straining  the  stresses  near  the  tip  of  an  infinitely  sharp  crack  are  (Refs.  3,4) 


ih 


-5- 


(18) 


-Km  sin  (0/ 2) 


K,,,  cos  (6/2) 


\j  2ttt 


°yz 


2nr 


or,  in  polar  coordinates 


K,„  sin  (0/2) 


K,,,  cos  (0/2) 


yj  2m 


(19) 


These  expressions  may  be  used  for  the  stresses  aV)  in  the  region  R2  -  The  corresponding  strain  energy 
density  is 


=  K2m/4Gitt 


(20) 


so  that  the  strain  energy  in  region  R2  is 


/  ^  CTjj  Cjj  dR  =  /  /"  (K^,/4G7rr)  r  d0  dr 

D  « 


o  -» 

=  K 6  /  2G 


(21) 


On  comparing  (21)  with  (9)  and  remembering  that  G  replaces  M  for  mode  III  straining,  we  find  that 

f,  =  %  (22) 


2.4.2  The  Factor  f2 

In  order  to  find  f2  the  strain  energy  of  the  stress  field  a-  must  be  calculated.  To  carry  out 
the  calculation  exactly  for  specimens  of  arbitrary  shape  would  be  difficult  if  not  impossible  and 
fortunately  is  not  necessary.  Upper  and  lower  bounds  on  the  strain  energy  can  be  established  by  using 
the  theorems  of  minimum  potential  energy  and  minimum  complementary  energy.  The  bounds  turn 
out  to  be  close  enough  to  fix  the  strain  energy  to  an  adequate  degree  of  approximation. 

Application  of  the  principle  of  minimum  potential  energy  requires  the  construction  of  a 
suitable  approximate  but  compatible  displacement  field,  while  application  of  the  principle  of  mini¬ 
mum  complementary  energy  requires  the  construction  of  an  approximate  stress  field  which,  however, 
must  satisfy  the  requirements  of  equilibrium.  Since  the  stress  field  is  concentrated  around  the  hole, 
suitable  approximate  fields  can  be  constructed  in  the  following  way.  With  centre  at  the  centre  of  the 
hole  draw  the  largest  possible  circle  which  lies  entirely  within  the  specimen  and  does  not  extend 
beyond  the  opposite  end  of  the  crack.  Examples  of  such  circles  are  shown  in  Figure  4.  Let  the  radius 
of  the  circle  be  c.  It  may  be  assumed  c  is  large  relative  to  the  hole  radius  b.  In  the  region  outside  the 
circle  of  radius  c  the  stress  and  displacement  will  be  taken  to  be  zero.  Thus  the  required  boundary 
conditions  on  the  exterior  of  the  specimen  are  satisfied.  The  stresses  and  displacements  inside  the 
circle  will  be  determined  so  as  to  satisfy  all  the  requirements  of  equilibrium  and  compatibility  as  well 
as  the  boundary  conditions  on  the  surface  of  the  hole  and  on  the  faces  of  the  crack.  If  in  addition  the 
stresses  satisfy  the  boundary  condition  of  zero  traction  on  the  circle  of  radius  c  then  we  will  have 
obtained  a  stress  field  which  satisfies  .the  equilibrium  conditions  throughout  the  entire  specimen.  By 
virtue  of  the  theorem  of  minimum  complementary  energy,  the  strain  energy  calculated  from  this 
equilibrium  stress  field  overestimates  the  correct  strain  energy.  On  the  other  hand,  if  the  displacement 
field  within  the  circle  is  made  to  satisfy  the  boundary  condition  of  zero  displacement  on  the  circum¬ 
ference  then  we  will  have  obtained  a  displacement  field  which  satisfies  the  requirements  of  compati¬ 
bility  throughout  the  entire  specimen.  By  virtue  of  the  theorem  of  minimum  potential  energy  the 
strain  energy  calculated  from  this  displacement  field  underestimates  the  correct  value. 
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ux  =  uy  =  0 


uz  =  u(x,y) 


It  follows  that  the  stresses  are  given  by 


oxz  =  G  9u/9x 


oyz  =  G  9u/9y 


®xx  ®yy  ®zz  -  ®xy  -  ® 


The  stresses  must  satisfy  the  equilibrium  equation 


from  which  it  follows 


9oxz 

-  + - =  o 

9x  9y 


V2  u  =  0 


At  this  point  it  is  convenient  to  change  to  polar  coordinates  because  they  are  more  sui‘ 
to  the  geometry  of  the  hole.  Equation  (26)  in  polar  coordinates  is 


92u  1  9u  1  92u  _ 

V2U  =  -  +  —  +  — r  — ;  —  0 

9r2  r  3r  r2  9 02 


while  the  formulas  for  stresses  in  polar  coordinates  are 
orz  *  G  9u/9r 


a0}  =  G  9u/(r90) 


The  boundary  traction  Tz  on  the  surface  of  the  hole  can  be  found  from  (19),  from  which  it  follows 
that  the  boundary  condition  for  the  stress  field  Oy  at  the  surface  of  the  hole  is 

K,u  sin  (0/2) 


y  2irb 


It  is  easily  verified  that 


u  =  (AVr  +  B/Vr)sin  (0/2)  (30) 

is  a  solution  of  (27)  for  any  values  of  the  constants  A  and  B,  and  further  this  solution  leaves  zero 
traction  on  the  faces  of  the  crack.  By  suitable  choice  of  the  constants  A  and  B  the  solution  can  also  be 
made  to  satisfy  the  boundary  condition  on  the  hole  and  the  chosen  boundary  condition  on  the  circle 
of  radius  c.  The  solution  which  gives  zero  traction  on  the  circle  of  radius  c  is  found  to  be 


u  *  \j2lrr  (K|U  6/G)  (1  +  r/c)  (1  -  b/c)*1  sin  (0/2) 


while  the  solution  which  gives  zero  displacement  on  the  circle  of  radius  c  is  found  to  be 


u  -  yj 2/xr  (K|||b/G)  (1  -  r/c)  (1  +  b/c)'»  sin  (0/2) 


The  strain  energy  calculated  from  the  former  solution  is 


The  strain  field  associated  with  mode  III  straining  is  of  the  type  known  as  “anti-plane” 
strain.  Anti-plane  strain  is  an  out-of-plane  shearing  action  whose  characteristics  are  that  the  only  non¬ 
zero  component  of  displacement  is  the  z-component,  and  that  this  component  is  independent  of  z. 
Thus,  in  die  usual  notation, 
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Km  b  (1  +  b/c) 
2G  (1  -  b/c) 

while  the  strain  energy  calculated  from  the  latter  solution  is 

KUI  b  (1  -  b/c) 
2G  (1  +  b/c) 


(33) 


(34) 


As  already  explained,  these  two  values  bracket  the  correct  strain  energy  of  the  field  o{ ..  Since  c  may  be 
assumed  large  relative  to  b,  it  follows  from  (33)  and  (34)  that  we  have  approximately 

/  \  olj  e'j  dR  =  K2m  bl 2G  (35) 

R.  2 

to  within  an  error  of  order  b/c  or  less.  On  comparing  (35)  with  (10)  we  conclude  that 


f2  =  >/2 


(36) 


2.5  Equivalent  Crack  Length  for  Mode  III  Straining 

As  we  have  already  seen,  terminating  a  crack  by  a  small  hole  increases  its  effective  length 
by  the  amount  2(f ,  +  f2 )  b.  In  view  of  the  values  obtained  for  f,  and  f2  this  amount  is  26  for  mode 
III  straining.  Since  the  crack  tip  in  the  ideal  specimen  was  taken  at  the  centre  of  the  hole,  the  tip  of 
the  effective  crack  of  the  real  specimen  therefore  lies  a  distance  b  beyond  the  edge  of  the  hole,  where 
b  is  the  radius  of  the  hole. 
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